Spectrum of anisotropic exponents in hydro dynamic systems with pressure 
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We discuss the scaling exponents characterizing the power-law behavior of the anisotropic com- 
ponents of correlation functions in turbulent systems with pressure. The anisotropic components 
are conveniently labeled by the angular momentum index I of the irreducible representation of the 
SO(3) symmetry group. Such exponents govern the rate of decay of anisotropy with decreasing 
scales. It is a fundamental question whether they ever increase as I increases, or they are bounded 
from above. The equations of motion in systems with pressure contain nonlocal integrals over all 
space. One could argue that the requirement of convergence of these integrals bounds the exponents 
from above. It is shown here on the basis of a solvable model (the "linear pressure model"), that 
this is not necessarily the case. The model introduced here is of a passive vector advection by a 
rapidly varying velocity field. The advected vector field is divergent free and the equation contains 
a pressure term that maintains this condition. The zero modes of the second-order correlation func- 
tion are found in all the sectors of the symmetry group. We show that the spectrum of scaling 
exponents can increase with i without bounds, while preserving finite integrals. The conclusion is 
that contributions from higher and higher anisotropic sectors can disappear faster and faster upon 
decreasing the scales also in systems with pressure. 
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I. INTRODUCTION 



Turbulent flows are often forced in an anisotropic fash- 
ion. The anisotropy has a significant effect on a variety 
of measured turbulence characteristics. We are interested 
in the effect of anisotropy on statistical quantities, espe- 
cially the low-order structure functions of velocity differ- 
ences across a scale R. In perfectly isotropic systems such 
objects are expected to display pure scaling behavior in 
the limit of high Reynolds number. We have suggested 
recently |0] that in the presence of anisotropy, the struc- 
ture functions are no longer pure power laws. Instead, 
components of the structure functions that belong to dif- 
ferent irreducible representations (sectors) of the SO (3) 
group possess different scaling exponents. Each of these 
sectors is characterized by the angular momentum indices 
£ and m. By projecting the structure function onto the 
different sectors, we could measure the universal 

scaling exponents in each sector separately. 

The spectrum of anisotropic exponents is particularly 
accessible in linear problems like passive scalar |^ and 
passively advected magnetic fields In both the mod- 
els, the equations of motion are isotropic and as a result 
the existence of universal anisotropic exponents can be 
proven Additionally the isotropy of the equations 

implies that the scaling exponents depend on £ but not 
on m. One of the important results of the analysis is that 
the discrete spectrum of anisotropic exponents is strictly 
increasing as a function of ^. This explains the isotropiza- 
tion of the statistics as smaller and smaller scales are 
observed. Since the scaling exponents C appear in power 
laws of the type (i?/L)'>, with L being some typical outer 
scale and R <^ L, the larger is the exponent, the faster is 
the decay of the contribution as the scale R diminishes. 
Therefore the gap between the leading, isotropic expo- 
nent, and the next available exponent governs the rate of 



isotropization. 

Experiments and simulations on Navier-Stokes (NS) 
turbulence also indicate that anisotropic sectors possess 
larger scaling exponents than the isotropic sector [^-Q. 
However, to date, the exponents for £ > 2 were not de- 
termined with sufficient accuracy. We thus do not know 
whether the higher sectors are characterized by an ever 
increasing exponents, or whether the exponents saturate. 
This issue is theoretically puzzling because of the effects 
of pressure. The inversion of the pressure in terms of the 
Green's function of the Laplacian operator introduces in- 
tegrals over the domain of turbulence. These integrals 
manifest the nonlocality of the problem, and are present 
in both the dynamical equation, and in the equations for 
the correlation functions. When considering a spatially 
homogeneous turbulence, the turbulent domain is usually 
taken to be infinite. The physical boundary conditions 
at scale L, are mimicked by employing a homogeneous 
forcing that acts at that scale. In this case, the integrals 
that result from the pressure term are over all M.^, and 
their convergence has to be guaranteed. The question 
that we want to address in this paper is the following: 
does the requirement of the convergence of the integrals 
necessarily bound the spectrum of the scaling exponents 
from above? Since the correlation functions appear in 
the integrand, an unbounded spectrum implies a rapidly 
increasing integrand as a function of the length scale. On 
the face of it, at some point the integrals must diverge in 
the infrared. It would appear therefore that either there 
must be a limit to the magnitude of the scaling expo- 
nents, or that the integrals converge due to an infrared 
crossover in the correlation functions. The latter sce- 
nario looks physically reasonable, yet in the presence of 
pressure integrals, seems to break scale invariance in the 
inertial range. To demonstrate that, consider a typical 
integral term of the form. 
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dyG{r-y)C{y) 



(1.1) 



with the same riddle. In Sec. IV we display the exact so- 



Here G{r) — — l/(47r|r|) is the infinite domain Green 
function of the Laplacian operator, and C(r) is some 
statistical object, which is expected to be scale invariant 
in the inertial range. If C(r) has an infrared crossover 
at scale L (or equivalently, the integral has an infrared 
cutoff at scale L), then the above expression will not be 
a pure power law of r, not even inside the inertial range. 
Then how is it possible that such an expression will can- 
cel out a local term of C(r), as is required by the typical 
equations of motion? 

This puzzle has led in the past to the introduction of 
the concept of "window of locality" Pp^]. The window 
of locality is the range for the scaling exponents in which 
no divergence occurs, even if the crossover length L is 
taken to infinity. For these exponents, integrals of type 
(1.1) are dominated by the range of integration y k r 
and are therefore termed "local" . In a "local" theory no 
infrared cutoff is called for. 

In this paper we show that scaling behavior of the cor- 
relation functions together with finite integrals over an 
infinite domain do not necessarily imply a bounded spec- 
trum of anisotropic exponents. Our strategy in this paper 
is to come up with a tractable example of a linear model 
with pressure, see Eq. (2.1). We refer to this model as the 
"linear pressure model". We approach the solutions of 
this model in two steps. First we distill yet another, sim- 
pler, exactly solvable model, which still poses the riddle 
of the Navier-Stokes problem. The exact solution reveals 
that the spectrum of scaling exponents is unbounded, 
and the convergence of the integrals is nevertheless not 
compromised. In the second step wc find the scaling ex- 
ponents of the "zero modes" 1lO| , |l l| of the linear pressure 
model, and use the conclusions of the simpler model to re- 
late them to the full solution. We show that also in this 
case the spectrum does not saturate in the anisotropic 
sectors. 

The linear pressure model and its simplified version re- 
veal two mechanisms that allow an unbounded spectrum 
of scaling exponents. First, a careful analysis of the win- 
dow of locality in the anisotropic sectors shows that it 
widens as £ increases. We always have a leading scaling 
exponent within the window of locality. Second, there is 
a more subtle mechanism that comes to play when sub- 
leading exponents exist outside the window of locality. In 
these cases we show that there exist counter-terms in the 
exact solution (not the zero modes!) that maintain the 
locality of the integrals. The bottom line is that in these 
models the anisotropic exponents are unbounded from 
above leading to a fast decay of the anisotropic contribu- 
tions in the inertial range. 

In Sec. we introduce the linear pressure model and 
derive the equations for the two-point correlation func- 
tion. We arrive at the form containing the dangerous 
integrals, and discuss again the fundamental riddle. In 
Sec. Ill we construct a simpler, exactly solvable model 



lution of the model, and discuss the windows of locality 
and the existence of an unbounded spectrum. In Sec. ^ 
we go back to the Linear Pressure model and offer a solu- 
tion of its zero modes. Section ^ offers a summary and 
a discussion. 



II. LINEAR PRESSURE MODEL 

A. The model 

The linear pressure model captures some of the aspects 
of the pressure term in Navier-Stokes turbulence, while 
being a linear and therefore, much simpler problem. The 
nonlinearity of the Navier-Stokes equation is replaced by 
an adverting field w{x,t) and an adverted field v{x,t). 
The adverting field w{x^t) is taken with known dynamics 
and statistics. Both the fields are assumed incompress- 
ible. The equation of motion for the vector field i;"(a;, t) 
is: 



dav" = 
c)„w" = 



(2.1) 
(2.2) 
(2.3) 



In this equation, /(a;, t) is a divergence free forcing term 
and K is the viscosity. The domain of the system is taken 
to be infinite. Following Kraichnan's model for passive 
scalar we choose the adverting field ■w{x,t) to be a 
white-noise Gaussian process with a correlation function 
that is given by 



S{t' - 1)0"!^ {r) = (u>"(a; + r, t')w^[x, t)) 
X"'3(r) = D"'3(r)-i:i"'3(0) 



(e + 2)^"^'-^ 



(2.4) 
(2.5) 

(2.6) 



The forcing f{x,t) is also taken to be a Gaussian white 
noise process. Its correlation function is 

F^f'ir/mt - t') EE irix + r, t)f{x, t')) . (2.7) 

The forcing is responsible for injecting energy and 
anisotropy to the system at an outer scale L. We 
choose the tensor function F"^{x) to be analytic in x, 
anisotropic, and vanishing rapidly for \x\ ^ 1. Analytic- 
ity is an important requirement. It means that F°'^{x) 
can be expanded for small \x\ as a power series in x"; 
as a result its leading contribution in the H. sector is pro- 
portional to given by d°'d'^x^Yhn{x). To see that 
this is the leading contribution the reader can consult the 
general discussion of the construction of the irreducible 
representations in Ref. All other analytic contribu- 
tions contain less derivatives and are therefore of higher 
order in x. 

'In order to derive the statistical equations of the corre- 
lation function of v°'{x, t), we need a version of Eq. (2.1) 



2 



without the pressure term. Following the standard treat- 
ment of the pressure term i n N avier-Stokes equation, we 
take the divergence of Eq. (2.1) and arrive at 



(2.8) 



The Laplace equation is now inverted using the Green 
function of infinite domain with zero-at-infinity bound- 
ary conditions, 

p{x) ^- JdyG{x-y)d,d^w''{y)v''{y) , (2.9) 



AttIxI 



(2.10) 



With this expression for p{x), Eq. (2.1) can be rewritten 
as 

dtv°'{x,t)+w>'{x,t)d^,v°'{x,t) (2.11) 
JdyGix - y)d,d^w^'{y)v%y) (2.12) 
-Kd'^v°'{x,t) = r(x,t) . 



B. Equations for the second-order correlation 
function 

The statistical object that we are interested in is the 
two-point correlation function of the field v"{x), 



C"^(r) EE {v"{x + r)v^{x)) 



(2.13) 



The last equation is identical to the equation for the 
second-order correlation function in the usual Navier- 
Stokes turbu lence, provided that w^^ is replaced with 
in Eq. ( 2.16 ). Indeed, the vexing problem that we face is 
being made very clear: if the triple correlation function 
has a power law dependence on r with an arbitrarily large 
exponent, how can the integral converge in the infrared? 
One possibility is that the scaling exponent of T"^{r) 
is sufficiently low, making the integral convergent. The 
other possibility is that the correlation function is scale 
invariant only in the inertial range and vanishes quickly 
after that, which is equivalent to the introduction of an 
infrared cutoff. However, the integral terms in the equa- 
tion probe the correlation function throughout the entire 
space. Therefore, a crossover behavior of the correlation 
function at the outer scale L, seems to contradict a pure 
scaling behavior of the correlation function in the iner- 
tial range itself. This in turn implies the saturation of 
the anisotropic scaling exponents. 

To proceed, we use the fact that the field w{x,t), as 
well as the forcing, are Gaussian white noises. This en- 
ables us to express T°'^{r) and the correlation of the force 
in terms of C°'^{r) and F°'^{r). In appendix ^ we use 
the well-known method of Gaussian integration by parts 
]l3[ that leads to the final equations. 



We find its equation of motion in two steps. First, we 
take the time derivative of C^^lr) using Eq. ( ^.ll| ), 

dt {v'^ix + r)v'^{x)) + {v^'ix + r)w''{x)d^,v'^{x)) (2.14) 
+ {v'^{x)w>'{x + r)df,v°'{x + r)) (2.15) 

- (v'^ix + r)af^) JdyG{x - y)d^d,w^'{y)v'' {y) ^ 

vf\x)df^^,~^ JdyG{x + r^ y)d^d,w'^{y)v''{y) ^ 

-K (v^'ix + r)d^v^\x)) - K {v^{x)d^v°'{x + r)) 
= {v^ix + r)f{x)) + {v^\x)r{x + r)) . 

To simplify the equations we define an auxiliary function 

T°'l^{r) EE a^*-) {v"{x + r)w^'{x)v^\x)) . (2.16) 

Using this definition and the space homogeneity of the 
statistics, we arrive after some algebraic manipulation to 
the following equation: 

5tC"''(r) - T-^f^r) - T^°'{-r) (2.17) 
+ jdyG{r-y)d^d,T'-''{y)+ jdyG{-r - y)d''d,T^%y) 
-2nd'^G"'^{r) = {v"{x + r)fl^{x)) + {v^\x)f°'{x + r)) . 



- JdyGir-y)d^d,T''''iy) 

~ jdyG{^r-y)d''d,T^'-{y) 
+2nd^G'^'^{r)+F"'^{r) , 



(2.18) 



(2.19) 



dZ. dyG{r - y)dr Ik^" iy)d^d,G-^ (y) 



1 

2' 
1 



-- \dyG{y)df'dr K>^''{y)d^d,C^^{r~y) 



These equations have to be supplemented with two more 
equations that follow directly from the definition of 
C"^(r), 



a„C"^(r) =0 , 

C"^(r) = C^"(-r) 



(2.20) 
(2.21) 



The first equation follows from the incompressibility con- 
straint of the vector field v{x^t), while the latter fol- 
lows from space hom ogeneity. Finally, we note that 
Eqs. ( 2.18| ) and ( 2.19| ) can be interpreted in a transpar- 
ent way, utilizing two projection op erato rs that maintain 
the right-hand side (RHS) of Eq. (2.18) divergence-free 
in both indices. To define them, let us consider a ten- 
sor field X"-^{r) that vanishes sufficiently fast at infinity. 
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Then the two projection operators and "Pr are defined 
by 

Vi^X'^^ir) EE X"^(r) - j dyG{r - y)a^X^^(y) , (2. 

Vj^X'^^ir) = X"^(r) - y dy G(r - y)d^X'^^{y) . (2. 

We observe that Vi^X°'^ and V^X"^ are divergence-free 
in the left and right indices, respective ly. Using these 
operators we can rewrite Eqs.( ^!lq - |2.19 ) in the form 



(2.24) 



T°'P{r) 



(2.25) 



dyGiy)^^^r\K^''{y)^^^,.C''^r - y) 



The projection in Eq. ( ^.25 ) guarantees that T°'^{r) is 
diver genc e free in its left index, while the projection in 
Eq. ( ^.24| ) guarantees divergence freedom in the right in- 
dex. 

Not all the terms in these equations are of the same 
nature. The integrals due to the projection operator are 
easy to deal with by applying a Laplacian on them. For 
example. 



(2.26) 



On the other hand, there seems to be no way to eliminate 
the last integral in Eq. (2.25), and therefore we shall re- 
fer to it as the "nontrivial integral" . Only for ^ = and 
^ = 2 it trivializes: the integral vanishes when ^ = and 
is proportional to C"^{r) when £, = 2. Unfortunately, in 
these extreme cases also the projection operator trivial- 
izes, and the effect of the pressure cannot be adequately 
assessed. We prefer to study the problem for a generic 
value ^ for which the incompressibility constraint and the 
pressure terms are nontrivial. 

We deal with the this problem head-on in Sec. 0. Due 
to the nontrivial integral, we will not be able to provide 
a full solution of C"'^(r), but only of the zero modes. 
However, before doing so we would like to study a model 
that affords an exact solution in order to understand in 
detail the issues at hand. In the next section we therefore 
consider a simplified model of the linear pressure model, 
yet posing much of the same riddle. 



III. AN EXACTLY SOLVABLE TOY MODEL 



We cons truct a toy model that inspired by Eqs. ( 2.1g ) 
and, ( 2.ig| ) for the correlation function in the linear pres- 
sure model. Within this model we demonstrate the strat- 
egy of dealing with the nonlocal-pressure term. Since it 
is a simplification of the statistical equation of the linear 
pressure model, the toy model has no obvious underlying 
dynamical equation. 



A. Definition of the toy model 

22) In the toy model, we are looking for a "correlation 
function" C"(r), whose equations of motion are 

5tC"(r) - -if^"^(r)a^5,C"(r) - 8^^^ (3.1) 

X JdxG{r - x)drKf'"{x)df,d^C^{x) 

+ Kd^C°'{r)+F"{r/L) , 
a„C"(r)=0. (3.2) 

Here F°'{x) is a one-index analog of the correlation func- 
tion of the original forces F°'^{x). Accordingly, we take 
it anisotropic, analytic in and rapidly vanishing for 
a;| ^ 1. As in the previous model, also here analytic- 
ity requires that the leading contribution for small |a;| is 
proportional to d"' x'^Ytrn{x) in the I sector. Accordingly 
it is of order x^^^ . 

The toy model is simpler than the Linear Pressure 
model in two aspects: First, the "correlation function", 
C"(r) has one index instead of two and therefore can be 
represented by a smaller number of scalar functions. Sec- 
ond, the unpleasant nontrivial term of the linear pressure 
model is absent. This will allow us to solve the model ex- 
actly for every value of ^. Nevertheless, the toy model 
confronts us with the same conceptual problems that ex- 
ist in the linear pressure model and in NS turbulence: 
can a scale invariant solution in the inertial range with a 
crossover to a decaying solution at scale L, be consistent 
with the integral term? If not, is there a saturation of 
the anisotropic exponents? 



Equation 3.1 can be rewritten in terms of a new pro- 



jection operator V, which projects a vector X"(r) on its 
divergence free part 



+ Kd^C"" + F° 



(3.3) 



where 



PX"(r) = X^{r) - 9" dy G{r - y)d^X''{y) . (3.4) 



We shall solve this integro-differential equation by first 
turning it into a partial differentail equation (PDE) using 
the Laplacian operator, and then turning it into a set of 
decoupled ODE's using the S0(3) decomposition. 

As in the linear pressure model, the nonlocality of the 
projection operator can be removed by considering a dif- 
ferential version of it 



d^VT°'{r) = d^T"{r) - a"9^r^(r) 



(3.5) 



In stationary condition 9tC" — 0, and therefore the dif- 
ferential form of the toy model is given by 
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= d^K''''{r)df,d^C°'ir 
a„C"(r) = . 



d°'drKt'''{r)d^,d^C^{r) 



(3.6) 
(3.7) 



We have reached a hnear PDE of order four. This PDE 
wiU be solved by exploiting its symmetries, i.e., isotropy 
and parity conservation, as demonstrated in the next sub- 
section. 



B. The SO(3) decomposition 



Equation 3.6 and the incompressibility condition of 
C"(r) are both isotropic and parity conserving. There- 
fore, if we expand C"(r) in terms of spherical vectors 
with a definite behavior under rotations and under re- 
flections, we would get a set of decoupled ODE's for their 
coefficients. For each sector {i,m), ^ > of S0(3) we 
have three spherical vectors: 



A?(f) 



(r) 
(r) 



(r) 



(3.8) 



Here ^^^(r) — r^Yim{i'), and see 0] for further details. 
The first two spherical vectors have a different parity 
than the third vector, hence the equations for their coef- 
ficients are decoupled from the equation for the third co- 
efficient. In the following, we shall consider the equations 
for the first two coefficients only, as they have a richer 
structure and larger resemblance to the linear pressure 
model. Finally, note that the isotropic sector, i.e., ^ = 0, 
is identically zero. To see why, notice that in this special 
sector there is only one spherical vector, A"(f) = r~^r'^ . 
Hence the isotropic part of C"(r) is given by c{r)r~'^r°' , 
c{r) being some scalar fu nctio n of r. But then the in- 
compressibility condition (3.2) implies that c(r) ~ r^"^ , 
which has a UV divergence. We therefore conclude that 
c(r) = 0, and restrict our calculation to ^ > 0. 

By expanding C"(r) in terms of the spherical vectors 
Ai and we obtain a set of ODEs [decoupled in the 
(^, m) labels] for the scalar functions that are the coeffi- 
cients of these vectors in the expansion. The equations 
for these coefficients can thus be written in terms of ma- 
trices and column vectors. To simplify the calculations, 
we find the matrix forms of the Kraichnan operator and 
of the Laplacian of the projection operator separately, 
and only then combine the two results to one. 



1. The matrix form of the Kraichnan operator 

To obtain the matrix of the Kraichnan operator in the 
basis of Ai and A2, we expand C"(r): 



C"(r) = ci(r)A?(f)+C2(r)AJ(f) 



(3.9) 



Using the basic identities of the ^tmi'i') functions (see 

i), 

d^'^imir) = , 
r-^a^$^m(r-) =M,,^(r) , 



a short calculation yields 
/CC"(r) = A''"'(r)a^a,C"(r) 



(3.10) 



2c'/ + 2{2 + 0-^-{2 + Oi£ + + 2)^J A?(r) 
24' + 2(2 + 0^ + 2(2 + 0^ 

C2 



-2(2-fe)^(^-l)-±jA^(r) . 

Therefore, in matrix notation, the Kraichnan operator 
can be written as: 



IC 



2Dr^ 



+2D{2 + Cy 
-D{2 + 0r^ 
r^K2 









(9 








(4 



_2 f {e + m + 2) 





£{l - 1) 



+ r 



(3.11) 



2. The matrix form of the Laplacian of the Projection 
operator 



Let 

T"(r)=ti(r)A?(f)+t2(r)A^(r) 
and applying a Laplacian to VT'^ , we get 



r r 



(3.12) 



(3.13) 



4'-^ + (2-^)^ 

Hence in matrix notation, 
ti\ f ~£ 



^2 



t" 

_ 1 

r\-\ 2-1 )[t'2 

i_ [£{£ + !) £{e-i){e + i) 

r2 y 
1 



1 



(3.14) 
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C. The matrix form of the toy model 

Now that the matrix forms of the Kraichnan opera- 
tor and of the Laplacian of the projection operator have 
been found, we can combine these two results to find the 



matrix form of the left-hand side (LHS) of Eq. (3^). To 
this aim let us define 



and from Eqs. ( [3.11 ) and (3.14) we get 



-1 
„(2) 



,(4) 
'1 

,(4) 



„(3) 

-1 

„(3) 



(3.15) 



(3.16) 



be omitted, thus reaching the following equation for the 
column vector (ci,C2): 




Finally, also the incompressibility constraint 
daC°'{r) = 0, can be expressed as a relation between 
ci(r) and C2(r), 



c'l + 2^ + ^4 - l{i - 1)^ = 



(3.21) 



This constraint has to be taken into account when solving 
Eq. ( p^ . 



Cl 
C2 



where the number in parentheses denotes the order of the 
derivative. The matrices Mi are given by 



^0 



J14 = Jr2i]^2 , 
I3 = 2eP2K2 

I2 = e(e - 1; 

Ii = (e-l)(^-2)P2Ki + 2(e-2 
+ (e-l)PiIKi-hPiKo + Po"" 
lo = (e - 2)(e - 3)P2Ko + (e - 2)Pi]Ko + PoKo 



(3.17) 



^2^2 + 2(e - l)P2l 
^Ki + P0K2 , 



2^0 



^1 



To find the RHS of Eq. (3.6) we expand the "forcing" 
F°'{r) in terms of the spherical vectors Ai and A2, 



F"(r) = /i(r)A?(f)-H/2(r)A^(f) 



(3.18) 



and applying a Laplacian we find the matrix form of 
(92F"(r), 



.f'{ + lf[-{i+m + '2)^h 

/2' + ^/2 + M-^(^-i);^/2 



Pi 

P2 



(3.19) 



At this point it is worthwhile to remember that the forc- 
ing term F°'{r/L) is assumed to be analytic. As a re- 
sult for r/L ^ 1, its leading contribution in the {£,m) 
sector is proportional to d"r^Y£m{r) ~ r^~^ . However 
d'^F°'{r/L) is also analytic, and must therefore also scale 
like r^~^ for small r, instead of scaling like r^~^ , which 
could be the naive dimensional guess. 

To proceed we restrict ourselves to finding the solu- 
tion in the inertial range and beyond. In these ranges 
the dissipative term Kd'^d^C°'{r) is negligible and can 



IV. SOLVING THE TOY MODEL 



A. The general solution 



The solution of Eq. (3jC) is somewhat tricky due to 
the additional constraint (3.21). Seemingly the two un- 
knowns Cl (r]and_C2 (?') are over determined by the three 



equations ( 3.20 



3.21 ) , yet this is not the case for the two 



equations (3.2C) are not independent (when considered 
as two scalar equations, resulting from the two-column 
vectorial equation). To see that this is the case and find 
the solution, it is advantageous to work in the new basis 



dl = Cl + ^C2 , 

d2 = -2ci + i{£ - 1)C2 



(4.1) 



In this basis the incompressibility constraint becomes 
very simple. 



d2 = rd\ 



(4.2) 



allowing us to express ^2 and its derivatives in terms of 
di. To do that in the framework of the matrix notation, 
we define the transformation matrix U 



U; 



1 £ 

-2 e{i-i) 

1 ^ £{£~l) ~£ 
2 



£{£ + 1) 



1 



so that, 



? )=U( 

£12/ \ C2 



(4.3) 
(4.4) 

(4.5) 



The equations of di{r) are the same as the equations for 
Ci(r), with the matrices replaced by 



(4.6) 
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and the sources p^ replaced by r^+Vsdf ^ + r^V^df^ + r^^W^df^ + r^" ^2^^ ^ 



P*A^v(pA , (4.7) +r«-3^i4^) + r^-'V,d, = ) , (4.9) 



d2 


= rd['^ , 






= r4^) - 






= rd(-^' - 


^24^) 


4'^ 


= r4^) - 


H34^) 


4^) 


= r4^) - 


m4'^ 



Their exphcit values are given by, 
V^ = D 
Vi = D 



P2 J \P2 

Notice that a divergence free forcing F"{r) will cause '^'lere are two dimensional vectors given by, 
pl{r) and P2{r) to be related to each other in the same 

way that di{r) and d2{r) are related to each other, i.e., ^/^ = p^^ ^ ^ (4 iq) 

P*2=r{piy. (4.8) /^x /Q 

y4 = N4 4 +N3 : 

Next, we perform the following replacements: \ / V 

^3 = N3 3 ^ + N2 J 

^1 = Ni ( } ) + No ( J 



= i^O 



We get an equation written entirely in terms of the func- 
tion di{r) and its derivatives, (4-11) 



2 

i6 + 6e 

V -nl 16 + 4^ 

V3 ^ 1 _4^2 _ 4^ + 32^ + 8 - + 

, -4^2 4£ ^ + 24 - _ + 

V2- iy\ _g^^2 _ - 4^ - 48 + 22^2 _ 2^2£2 _ 3^2^ ^ 



, _(^ + 2)(-6^ + e£2+e^ + 4£2 + 4^) 

''^1 " I (e + 2)(6e2 - e^£2 _ ^2^ _ ^^2 _ ;l8^ _ ^£ + ^4 _^ ^^£2 _^ + 10^) 



Vo = D 



^(e + 2)(e-4)(£-l)(^ + 2)(£ + l) 



Equation 49 is for a column vector, and can be re- 
garded as two scalar differential equations that we refer 
to as the "upper" and the "lower" . The upper ODE is of 
the fourth order, while the lower ODE is of fifth order. 
Not surprisingly, the lower equation is the first derivative 
of the upper equation, provided that F°'(r) is divergence 
free. Hence the two equations are dependent, and we re- 
strict our attention to the upper equation. To simplify 
it, we divide both sides by Dr^ , replace di(r) by tp{r) 
and define the RHS to be the function S{r), 



S{r) = D 



(4.12) 



After doing so, we reach the following equation: 



r + '^3— 1- 02-^5- + 01^3- + 0.0^ = b(r) 



(4.13) 



Its homogeneous solution is easily found once we substi- 
tute. 



(4.14) 
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The scaling exponents are the roots of the polynomial, 

F(C)=C(C-l)(C-2)(C-3) (4.15) 
+ a3C(C - 1)(C - 2) + a^CiC - 1) + aiC + oo 

and are found to be real and nondegenerate. Two of them 
are positive while the other two are negative, given in a 
decreasing order by 



Ci = 


1 

^2 ^ 






1 

^2 ^ 




C3 = 


1 

^2 ^ 




C4 = 


1 

~2 ^ 


1 

2^ 


where 







A{£,0 + 


VB{£,0 , 


(4.16) 


A{i,0- 


VBi£,0 , 




^{£,0- 






A{£,0 + 






-2^ + 5 


+ 4£ + 4£2 , 


(4.17) 



B(£, = -8r£ - + i6r + sr^-" + (4.18) 

-8^^^ - 8^£ - 32^ + 16 + 6A£ + 64^^ . 



In the limit ^ ^ the roots become: 

C2 = ^ - 1 , 



(4.19) 



C3 
C4 



Figure 1 displays the first few exponents as a function 
of We note that the spectrum has no sign of satu- 
ration as £ increases. Before we discuss the meaning of 
this observation we will make sure that these solutions 
are physically relevant and participate in the full (exact) 
solution including boundary conditions. 

The general solution of Eq. ( 4.13| ) is traditionally given 
as the sum of a special solution of the non-homogeneous 
equation plus a linear combination of the zero modes. 
However, when attempting to match the solution to the 
boundary conditions, it is more convenient to represent 
it as 



4 

E 



-1 (0- Ci)---(C.-C4) 

^ V ' 

all different roots 



J rrii 



(4.20) 



where the free parameters of the solution are the four 
constants m^. Indeed a change in rrii is equivalent to 
adding to the solution a const x r^* . In the next subsec- 
tion we find the values of rrii that match the boundary 
conditions, and discuss the properties of the solution. 





— — ' T^^- 




£ = 6 - 




£ = 4' 




£ = 2 - 





42 



£ = 8 




£ = 6 




£ = 4' 




£ = 2 
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B. Boundary conditions and inertial-range behavior 



From Eq. (4.20) it is clear that the only values of rrii 
that guarantee that the solution remains finite as r ^ 
and that it decays as r ^ c« are mi = m2 = +00, 

TO3 = 7714 = 



ip (r) 



(Ci-C2)(Ci-C3)(Ci-C4)y. 

T-C2 



(C2 


-Cl)(C2-C3)(C2 


-C4) 








(C3 


-Cl)(C3-C2)(C3 


-C4) 




J-C4 





(C4-Cl)(C4-C2)(C4-C3)7o 



dxx^-'^' S{x) 
dxx^~^^ S{x) 
dxx^-'^'> S{x) 
dxx^''^^ S{x) . (4.21) 



To unders tand the asymptotics of this solution we find 
from Eq. ( 4.12| ) that for a; -C L, S{x) has a leading term 
that goes like x^~^~^, whereas for for x :s> L, S{x) de- 
cays rapidly. It is now straightforward to prove that for 
r <^ L, the C3 and ^4 terms scale like r^+^~^ ^ the (,2 term 



scales like r''^, and the term scales like r'^'^ for values 
of ^ for which Ci < ^ + 3 — ^ and like otherwise. In 

addition, it is easy to see that for r ^ L, ip{r) exhibits an 
algebraic decay: the Ci and (2 terms decay rapidly due 
to the decay of S{x) whereas the (3 and C4 terms decay 
algebraicaly like r^' , respectively. The asymptotics of the 
full solution are thus given by 



ip{r) 



r < L 
r-> L 



(4.22) 



The obvious conclusion is that there is no saturation in 
the anisotropic scaling exponents as t increases. The lack 
of contradiction with the existence of an integral over all 
space has two aspects. The main one i s sim ple and obvi- 
ous. The integro-differential equation (3.1) for C" has a 
differential version (|3.6|) . Solving the differential version, 
we are unaffected by any considerations of convergence 
of integrals and therefore the solution may contain expo- 
nents that increa se wi th I without limit. Nevertheless, 
the full solution (4.21) exhibits a crossover at L: it in- 



creases in the inertial range r <^ L and decays for r ^ L. 
Thus plugging it back to the integro-differential equation 
we are guaranteed that no divergence occurs. 

The question why the cross-over length L does not 
spoil the scale invariance in the intertial range still re- 
mains. The answer is found in d iffer ential form of the 
equation of motion, given by Eq. (3^). From this equa- 



tion we find that the integrand is a Green's function times 
a Laplacian of a tensor. By definition, such an integral lo- 
calizes, i.e., it is fully determined by the value of the ten- 
sor at the external vector r. In the language of Eq. (1.1) 

The second and less obvious aspect is that the window 
of locality widens up with £. This is due to the can- 
cellations in the angular integration of the anisotropic 



solutions that are due to the orthogonality of the Y^mif") 
and their generalizations B^^^{r). To demonstrate this 



consider again the simple integral (|l.l|), and assume that 
C{y) belongs to (^, 777) sector, i.e.. 

Civ) = aiy)YerrAy) ■ (4.23) 
For y 3> 7", we may expand the Green function in r/y, 
1 



G{r -y) = - 



47r|r — y\ 

1 °° 



(4.24) 



, ^ ■ y 

y 



Here a„ are Taylor coefficients. Obviously the dangerous 
terms for the infrared convergence are those with low val- 
ues of 77. However, all these terms will vanish for n < £ 
in the angular integration against Ypmiv)- The reason is 
that all these terms are of the form r"^y"2(^r ■ y)^^ with 
713 < £. The angular part here has projections only Y^'m' 
with £' < < £. The first term to contribute comes 
when n = £, and is proportional to the amplitude inte- 
gral J^dy y'^aim{y)y^^^^ ■ For a power law ai^iy) ^ y^ 
this implies locality for 



X<£-2 



(4.25) 



instead of A < —2, as in the isotropic sector. The lower 
bound of the window of locality is also extended and a 
similar analysis for y <C 7- leads to A > — — 3. For the 
toy model this translates to the window of locality 



£-^<Q<£ + l-^ 



(4.26) 



From the previous analysis we find that the leading power 
law of the full solution in the inertial range is r^^ , which 
is inside this "extended" window of locality. Neverthe- 
less, the subleading power r^^ originating from the first 
term in Eq. (4.21) is above this window, and its pres- 
ence in the solution can be explained only using the first 
mechanism. 

We will see when we turn back to the linear pressure 
model that both these mechanisms operate there as well, 
leading again to a lack of saturation in the exponents. 



V. SOLVING THE LINEAR PRESSURE MODEL 

We now return to t he li near pressure model. As men- 
tioned before in Sec. [IB, we see no way of eliminating 
all integrals from the equation and therefore we will not 
look for a full solution. Nevertheless, we shall be able 
to calculate the zero modes and hence the scaling ex- 
ponents. Our strategy relies heavily on the conclusions 
of the last section: we will apply two Laplacians to the 
equation for C"^(r) in order to eliminate the integrals of 
the two projection operators "Pl and 'Pr,. The resulting 
equation will still contain the nontrivial integral. Using 
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numerical integration we will solve the homogeneous part 
of this equation, i.e., we shall find its zero modes. These 
are scale invariant solutions that solve an equation con- 
taining an integral. Their exponent must therefore lie 
within the "extended" {i dependent) window of locality. 
Finally, we will argue that these zero modes are a part 
of the full solution that decays for r L, and therefore 
solve the original equation as well. 



A. Equations for the zero modes 



We start from Eqs. ( 2.24 ) and ( 2.25| ). In appendix ^ 
we perform integration by parts and algebrai c ma nipu- 
lations to bring the nontrivial integral in Eq. ( 2.25 ) to a 
more tractable form. The result of this process is 



1 12^1) 



2(^-3)(e-5) 



JdyG{y)i/-^d'C'^^{r-y) , (5.1) 



which is true for every ^ 1. The ^ = 1 case will not 
be treated here explicitly. Nevertheless, we will argue a 
posteriori that the results for ^ = 1 can be deduced from 
the C 7^ 1 results by continuity. 

Looking at Eq. (5.1), we note that when ^ = 2, the 
integral on the RHS of the above equation trivializes to 
a local term C"'^(r). In this limiting case the model can 
be fully solved utilizing the same machinery used in the 
previous section. The solution can then be used to check 
the zero modes computed below for arbitrary value s of ^ . 

To proceed, we substitute Eq. ( ^ ) into Eq. ( 2.24 ), 
noting that the projector leaves the nontrivial inte- 
gral in Eq. ( ^.1| ) invariant since it is divergence-free in 
both indices. Setting dtC°'^{r,t) = in the stationary 
case, we arrive at the following equation: 







PnAK'^^d^d^C'''' (r) 



(e-3)(e-5) 



jdyG{y)y^-^d^C--^{r-y) 



(5.2) 



As in the toy model, we apply two Laplacians to the 
above equation in order to get rid of the integrals of the 
projection operators and obtain 



a/3 



(r) 



12CD 



(e-3)(c-5) 



JdyGiy)yi-^d'C''^{r-y) 



(5.3) 



Here and in the following, the operator 9^" should be 
interpreted as (9^)". We now seek the homogeneous sta- 
tionary solutions of C"''(r) in the inertial range (zero 



modes). These satisfy the equations obtained by ne- 
glecting the dissipation and setting the forcing and time 
derivative to zero, 

^ (9^if^''(9^a^C"^(r) + d°'d'^drdaK''''d^d^C^''{r) (5.4) 



+ 



Let us now define the RHS of the above equation as the 
"zero-modes operator" 0{^) and write the zero- modes 
equation compactly as 



= 



d{^)C°'^] (r) . 



(5.5) 



To solve it, we write the solution C"^{r) in a basis that 
diagonalizes O(^). This is done in the next subsection. 

B. The SO(3) decomposition 

To diagonalize O(^) we must look for its symmetries by 
looki ng for the operations that commute with it. From 
Eq. ( |5.4| ) it is easy to see that these are rotations, scal- 
ing, permutation of indices, and flipping of r. As a result 
0{£_) is block diagonalized by tensors C°'^{r) that have 
the following properties: 

• They belong to a definite sector {i, m) of the SO (3) 
group. 

• They have a definite scaling behavior, i.e., are pro- 
portional to r^ with some scaling exponent A. 

• They are either symmetric or antisymmetric under 
permutations of indices. 

• They are either even or odd in r. 

In 1^ we discuss these types of tensors in detail. Here 
we only quote the final results. In every sector (^, m) 
of the rotation group with t > \, one can find nine in- 
dependent tensors X"^{r) that scale like r'^. They are 
given by r^B'^^ ^(f) where the index q runs from 1 to 9 
enumerating the different spherical tensors. These nine 
tensors can be further subdivided into four subsets: 

• Subset I of symmetric tensors with (— )^ parity, 
containing 4 tensors. 

• Subset II of symmetric tensors with {~Y~^^ parity, 
containing 2 tensors. 

• Subset III of antisymmetric tensors with (— )^'^^ 
parity, containing 2 tensors. 

• Subset IV of symmetric tensors with (— )^ parity, 
containing 1 tensor. 
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Due to the diagonalization of O(^) by these subsets, the 
equation for the zero modes fohates and we can compute 
the zero modes in each subset separately. In this paper, 
we choose to focus on subset I, which has the richest 
structure. The four tensors in this subset are given by 



(5.6) 



B2/m(r) = r-^[r"9^ + r/'a"]<i>,,„(r) , 

We expect the calculation of the other subsets to be eas- 
ier once one is famihar with the techniques we are about 
to develop. FinaUy. we note that since the correlation 
C°'^{r) has to fulfill C"^(r) = C^"{-r), our subset 1 
solution will be valid only for even is. 
Expanding C"^{r) in subset I, 



+ C3S3L(^)+C4i?4"L(^) 



Q/3 



(5.7) 



and plugging it back into PDE (5^), we obtain a linear 
equation for the coefficients ci,C2, 03,04 that depend on 
the parameters A, £, ^. In the four-dimensional space of 
column vectors (ci, C2, C3, C4), we can write it as 



O(A;^,0 





/Ci\ 


) 






C3 







, 



(5.8) 



where 0(A; ^) is a 4 x 4 matrix that represents the zero- 
modes operator 0{£_). 

To conti nue, w e note that due to the incompressibility 
constraint ( ^.2C| ) of C"'^(r), not all vectors (ci, C2, C3, C4) 
are allowed; for a given A, £, only certain combinations 
of the B"^^{r) lead to a divergence free C"'^(r). A 
simple calculation ||l| shows that these belong to a two- 
dimensional subspace, which is spanned by the "incom- 
pressible vectors" 



^i(A;^)) = 




(5.9) 



The zero modes exponents A can be found by requiring 
the equation 



^,0 fli "i(A;^)) +02 W2(A;^)) =0, (5.10) 



to admit a nontrivial solution. The explicit computa- 
tion of the matrix 0(A;£,^) is technically very cumber- 
some due to the presence of the integral term. We shall 
there fore use an implicit method to determine whether 
Eq. (5.10) has a nontrivial solution or not. The basic 
idea is t hat the calculation of the nontrivial integral in 
Eq. (5.4) can be simplified if we contract its free indices 
with the two isotropic tensors farp and Saf3- Therefore, 
instead of solving Eq. ( |S.4[ ) explicitly, we will contract 
its RHS with these two tensors and require that the two 
resultant scalars vanish simultaneously. Obviously this 
would provide us with a necessary condition for the solv- 
ability of Eq. (5.4). Nevertheless, we shall see that it is 
also a sufficient condition. 

Let us write the two tensors faf^ and Saf3 in matrix 

notation as two row vectors ( wi (i) 



W2{1) 



given by 



(1 2i 1 l{t — 1)) contraction with rafp , 

( 1 2f 3 ) contraction with 5ap ■ 

(5.11) 



The contraction of da[3 and f^f^ with another tensor is 
translated to the usual matrix multiplication of these 
row vectors with a column vector. For example, if 

cj is a column vector whose components are given by 

Ci, i — 1, ... 4, then 

Sc.(3{c,B'^Uf) + C2<™(^-) + c.B^Lir) + c,B2Ur)) 
= fci -f 2^C2 + 3c3)Y(,„(f) 

'w2{i) c)Yt^{r) . (5.12) 



Returning to the zero-mode equations, we now define 
the 2 X 2 "reduced matrix" Oij{\]l,i) by 



(5.13) 



Obviously, the existence of a nontrivial solution that 
makes the two contractions vanish is equivalent to the 
requirement that Oij{X;£,$,) is singular, i.e., 



(5.14) 



The above condition is also sufficient for the solvabil- 
ity of the zero-modes equ ation. To see that, notice 
that the RHS of Eq. (|5j) and therefore the RHS of 
Eq. (5^) produce tensors that are divergence- free in both 

indices. Thus the vectors 0(A;^,^) Ui{X;£)'^ will be- 
long to the two-dimensional subspace that is spanned by 
Uj (A -f ^ — 6; £) ). Since the transformation matrix 



(5.15) 



= (^w,{£)\uj{X + ^-6;£)) , 
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is nonsingular for all values of A and ^, we find that 
0(A;£,C) Mi(A;£)^ and 0(A; ^, U2(A; are linearly 
dependent if and only if Eq. holds. 

Looking at Eq. (^), we recognize that 0{£,) is a sum 
of a differential and an integral operator. Consequently, 
0(A; £, £,) and Oij{X; £, £,) can be written as a sum of cor- 



responding parts, 

0(A; e, C) = 0*-^(A; i, + 0*"*(A; i, C) , 
O.,(A;£,0 = O*^(A;f,a + O^f(A;^,0 ■ 



(5.16) 



These parts are of a different nature and thus we dedicate 
two different subsections for their calculation. 



1. Form of the differential operator 



The calculation of O*^ (A; £, 

can be done directly by calculating the matrix 0'^*-''(A; ^, employing the same 



techniques t hat a re used in Refs ||l|,^ and in sec. [II. Here we merely give the mid and final results. 
From Eq. (5.4), we see that the differential part of the zero-mode operator is given by 

To find its matrix representation it is convenient to first calculate following operators: 
• L(A;^) : X"/' — * d^X^'f^ 



(5.17) 



L(A;^) 



/ (A-^-2)(A+j + 3) 
2 (A-^)(A + £+l) 



2 

Vo 



ix-e){x + e + i) 

{X-£ + 2){X + i-l) / 



(5.18) 



l(A;^) : X'^P — > d°'d'^drd„X^° 



l(A;^) 



(A-£-2)(A-^-4) 
A-^-2 
A-^-2 



(A + l)(A + 2) 
2^(A + 2)(A-^) 
(A + l+^)(A-£) 
^(A-^)(^-l)(A-^-f2) 



T 



(5.19) 



• M2 



(A; I) : X"l^ — > d"drX^(^ + d^drX""^ 



2{X- e-2){X + 2) 2i{X - £){X - i - 2) 2{X ~ £ ~ 2){X ~ £) 

A + 2 (A-^)(^ + A + 3) 2(A-^) (A - ^)(^ - 1)(A - £ + 2) 

2(A + 2) 2^(A-^) 2(A-^) 

2(A + 3) 2 2{i - 1){X - £ + 2) 



(5.20) 



IK(A;^,e) : — > Kt"'{^)d^,d^X°''^ 



K(A; £,0^D {i + 2)L(A; £) - £,X{X - 1)1 



(5.21) 



With these four matrices, the matrix form of the dif- 
ferential part is written compactly as 



])*/(A;£,^)= L(A + e-4;£)L(A + e-2;£) (5.22) 
+ Mi(A-t-e-2;^) 
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Of*^(A;i?,^) is then computed directly from definition. 



2. Form of the integral operator 



The integral part of the zero-modes operator is given 



by 

[0"*(C)C"''](r) 

(e-3)(e-5) 



jdyG{r~y)\r-y\i'^d'C'-f\v) . (5.23) 



To calculate the reduced matrix ( ^.13 ), we have to com- 
pute the contraction of with f^f^ and Sap- 
To do that, le t th e tensor X°-^{r) be of the form of the 
trial solution (5.7), i.e., let it belong to the {i,m) sector 
of SO (3), be divergence free in both indices, and propor- 
tional to r^. If we denote in matrix notation X°'^{r) 

X ) , then from the isotropy of the 



by the column vector 
tensors fafp, Sap, we have the following useful identities: 

fafpX^^ir) = (wi{l)\xy^Y,^{r) , (5.24) 

SafiX^^ir) = (w2{e)\x)r^YUr) ■ 

Consider now the contraction of 

JdyG{r - y)\r - y\^-^X''^\y) , (5.25) 

with Sap 

h^Sap jdyG{r-y)\r-y\^-^X'^P{y) (5.26) 
= (w2{t)\x) JdyGir-y)\r-y\i-^y^Yi,n{y) ■ 

If A is in the window of locality of the integral above, 
then the integral will converge and be proportional to 
r^'^^Yim{r)- The scaling exponent A -f ^ follows from 
power counting [remember that G{r) ^ l/r] while the 
angular dependence is a result of the isotropy of the in- 
tegration over all space. This leads us to define the pro- 
portionality factor A{\;(^,£), 

A{X;£,Or^^^Yera{r) (5.27) 
= JdyG{r -y)\r- y\^-^y^Y,^{y) , 

with which we can write 

h = (w2{£) x)A{X;£,Or^+iYe,n{r) . (5.28) 



The prefactor A{X; £, ^) is the only part of the calcula- 
tion that cannot be done analytically. However, instead 



of calculating the integral in Eq. (5.27) numerically, we 
can expand the integrand as a Taylor series and write 
A{X;£,^) as an infinite series of poles in A. This is done 
in detail in Appendix ^ with the result that 



A{X;£,0 = 
1 °° 



(5.29) 



g=0 



1 



1 



A + 3 + ^ + 29 X + ^-£-2q 



The coefficients aq{£,^) are given in Appendix ^ Wc no- 
tice that the window of locality in the definition ( ^.27 ) of 
A{X] ^, £) can be identified from the positions of the poles 
in Eq. ( ^.29 ). Indeed, the boundaries of the window of 
locality are determined by the q — Q term and located at 
A = —£ — 3 and X = £ — When A hits these boundaries 



th e inte gral in Eq. ( 5.27 ) diverges corresponding to a pole 



5.29. However, the above formula is valid also for val- 



ues of A outside this window of locality for which the 
formal definition of A(A; £, ^) makes no sense. The rele- 
vance of these values to the full solution must therefore 
be addressed. This will be done in the next subsection 
where we present the values of the scaling exponents. 
Let us now consider the fafp contraction 



h = r-^ jdy G{r ^ y)\r - yf-^rarpX^\y) . 
Using the identity 

Tarp = {ra - ya){rp - yp) + {Va - ya)yp 

+ Vairp - yp) -f yaVp , 
we can decompose the integral into four terms: 
• The {ra - ya){rp - yp) term: 

Recalling that 

G(r - y) ^ 



(5.30) 



(5.31) 



47r|r — y\ 



(5.32) 



(5.33) 

?+2 



it is easy to verify that 

G{r -y)\r - y\^^'^{ra - ya){rp - yp) 
= j^^^^^^dadpGir - y)\r - y 

- -^z:\G(t -y)\r- y\^Sap ■ 



Plugging this into the integral, the term with the 
derivatives will vanish due to integration by parts 
and the divergence- free X°'^[r). We are therefore 
left with 

- y%G(r - y)\T - y\^apX'^^{y) (5.34) 



e-1 



W2{£) x)ldyG{r-y)\r-yh^Y,„,{y) 
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Using the identity 

d^y^+^YeUy) = {X + 2-£){X + 3 + i)y^Y,^{y) , 

(5.35) 

we further integrate by parts the last integral to 
finally obtain 



(A + 2-£)(A + 3 + ^) 



W2{t) x)A{\ + 2■^,S) 

(5.36) 



The {ra - Va)vp, Vairfs - Up) terms: 

Using the same tricks as in the previous term, 
i.e., integration by parts and the fact that X°'^{y) 
is divergence free, we can easily show that both 
these terms vanish. 



• The yaVp term: 



Using Eq. ( |5.24| ) and the definition of A(A;£,^), 
wc directly get 

iwi{i)\x)A{\ + 2-i,i)r^+^Yt^{r) . (5.37) 



Gathering all the terms, we see that the contraction with 
fafp is equal to 



Wl{l) 



(A + 2-£)(A + 3 + £) 



(5.38) 



To conclude, the matrix 0™*(A;£, ^) is given by the fol- 
lowing equations: 



u,(A;^,0 



12^D 



(c-3)(e-5) 



— L(A-4;^)L(A-2;^)L(A;^) u^{\■^)) , 
5) / 



(5.39) 



wi{e) w,(A;£,e) 



OZ' (A; e, = ^(A - 6; i, {w2 (i) \ ^.(A; £, o) 



(A-4-£)(A-3 + £) 



W2{£) u,{X;£,0 



C. Results 

Figure 2 shows the leading scaling exponents of the 
linear pressure model for £ = 0,2,4,6,8,10. The re- 



sults where obtained by numerically solving Eq. (5.14) 
for ^ = 0, 0.01, 0.02, . . . , 1.99, 2. The prefactor A{X; £, ^) 
was calculated using the formulas ( C12| ) and ( pl3| ) where 
the infinite series of poles was truncated typically after 
100 poles, when it was clear that relative contribution 
from the consecutive pole was smaller than the machine 
precision (about 10^^^). Additionally, the numerical re- 
sults were compared to the analytical results of ^ = 0, 2 
and of ^ = and £ — 2 and found to be correct within a 
relative error of lO^^*'. 

From Fig. 2, we see that in the isotropic sector and in 
the £ — 2 sector, the leading exponent is C = 0, corre- 
sponding to the trivial C"''(r) = const solution. These 
zero modes will not contribute to the second-order struc- 
ture function, which is given by 



(5.40) 



and so we have to consider the zero mode with the con- 
secutive exponent. In the isotropic sector this exponent 



is exactly ^ — 2 — ^ as can be proven by passing to Fourier 
space. This special solution is a finger-print of the exis- 
tence of a constant energy fiux in this model. Indeed 
just like in Navier-Stokes turbulence, one can show ana- 
lytically that the isotropic part of the triple correlation 

function (^v°' {x)w^ {x -t- r)v^{x + r)^ is proprtional to r 

and hence <S'^"f^g^(r) ^ r^^^. 

Returning to the main question of this paper, we see 
that no saturation of the anisotropic exponents occurs 
since the leading exponent in every £ > 2 sector is 
({t) ~ — 2. These exponents are within the window of 
locality of Eq. (|^, which is given by-^-3<C<^-C- 
However, the next-to-leading exponents (that are the 
leading ones in the structure function for £ = 0,2) are 
already out of this window and their relevance has to be 
discussed. We propose that the same mechanism that 
works in the toy model (see sec. [V B| ) also operates 
here and that all these higher exponents can be found 
in the full solution. To understand this, let us write a 
model equ ation for the correlation function in the spirit 
ofEq.(yl), 



VC{r) 



dvK{r-y)C{y) = F{r) 



(5.41) 
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with K being some kernel, and V being some local dif- 
ferential operator. In view of Eq. (5^), the differential 
operator T> should be regarded as the Kraichnan opera- 
tor and the integral term should be taken for all integral 
terms in the equation including integrals due to the pro- 
jection operators. These integrals create a window of 
locality that we denote by Ajq^ < A < Ajjj. Any pure 



scaling solution C(r) 



with A outside the window 



of locality will diverge and h ence will not solve the ho- 
mogeneous part of Eq. (5.41). Nevertheless, we will now 
demonstrate how this zero mode can be a part of a full 
solution without breaking scale invariance. For this we 



act with a Laplacian on both sides of Eq. (5.41) in order 
to get rid of the projection operators integrals. Of course, 
like in the linear pressure model, this will not eliminate 
all integral terms, and thus we can write the resultant 
equation as 



dyK{r - y)d^C{y) = d^F{r) . (5.42) 



Our main assumption, which was proven analytically in 
the simple case of the toy model, is that the above equa- 
tion has a solution that is finite for all r and deca ys fo r 
r ^ L. Let us now consider the zero modes of Eq. ( 5.42 ); 
their exponents have to be within the "shifted" window 
of locahty Aj^^ -I- 2 < A < Aj^j + 2. Suppose now that 
with Ajjj < A < Aj^j -I- 2 is such a sol ution , which is 
therefore part of the full solution of Eq. ( 5.42| ) . We now 
claim that this solution also solves the original equation 
[Eq. (5.41)], hence allowing the existence of scaling expo- 
nents outside its window of locality. To see that, we first 
notice that since the ful l solu tion decays for r 3> i, then 
all the integrals in Eq. ( ^.41 ) converge and are therefore 
well defined. All that is left to show is that the equation is 



indeed solved by C(r). But this is a trivial consequence 
of the uniqueness of the solution for Laplace equation 
with zero at infinity boundary conditions. Indeed, if we 
denote the integral term in Eq. (5.41) by 



/(r) = J dyK{r-y)C{y) , (5.43) 



then from Eq. ( 5.42| ) we have 

d^Iir) = d^[F{r)-VCir)] , 



(5.44) 



and since both /(r) and F{r) — 'DC{r) decay as r ^ oo, 
then they must be equal. Of course no breaking of scale 
invariance occurs because the equation is satisfied and 
F{r) — 'DC{r) is a sum of an inhomogeneous solution 
and power laws. 

Returning to the linear pressure model, we have shown 
that not only the first, leading exponents in every sector 
are legitimate but also the next few exponents. These 
exponents are inside t he s hifted window of locality of the 
"Laplaced" equation (5.4), which is given by — ^ -f 1 < 

At this point, we may ask whether this is also the case 
for the other exponents, which are outside this shifted 
window of locality. In light of the above discussion, it is 
clear that all of them may also be part of the full solution 
for we can always differentiate Eq. (5.2) sufficient number 
of times, thus shifting the window of locality to include 
any of these exponents. However, this procedure is un- 
necessary once we have written the prefactor A(\\ £, ^) as 
an infinite sum of poles in A. In that case the equation 
is defined for all values of A except for a discrete set of 
poles, enabling us to look for exponents as high as we 
wish. 



15 





7 
6 

5 

3 

1 


1 





0.5 



1.5 



5 

4 - 
3 
2 



0.5 



1.5 




10) 




FIG. 2. Leading scaling exponents for the first few Ps in the Linear Pressure model. The dashed line indicates the upper 
bound of the window of locality 



VI. SUMMARY AND CONCLUSIONS 

1 

The main question raised and answered in this paper 

is whether the existence of the pressure terms necessarily 
leads to a saturation of the scahng exponents associated 
with the anisotropic sectors. Such terms involve integrals 



over ah space, and seem to rule out the existence of an 
unbounded spectrum. We have discussed a mechanism 
that allows an unbounded spectrum without spoiling the 
convergence of the pressure integrals. The mechanism 
is demonstrated fully in the context of the simple toy 
model and we proposed that it also operates in the case 
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of the linear pressure model. The mechanism is based 
on two fundamental observations. The first one is that 
the window of locality widens up linearly in £ due to the 
angular integration. The second, and more important, is 
that a scaling solution with an unbounded spectrum can 
exist as a part of a full solution, which decays at infinity. 
Indeed pure scaling solutions cannot solve by themselves 
the zero-modes equation if their scaling exponent is out 
of the window of locality. However, the zero-modes are 
always part of the full solution that decays to zero once 
r ^ L and we have shown that if such a solution solves 
a differential version of the full equation, it must also 
solve the original equation. Therefore by differentiating 
the full equation sufficiently many times, we can always 
reach a differential equation with a window of locality as 
high as we wish. In that equation we can find zero-mode 
solutions with arbitrarily high exponents (notice that in 
the toy model, it was sufRcient to differentiate once to 
get rid of all integrals, thus obtaining an "infinitely wide" 
window of locality). But since these zero-modes are part 
of a full solution that decays at infinity, then this solu- 
tion is also valid for the original equation, hence showing 
that in the full solution there can be power laws with 
arbitrarily high exponents. 

Finally, we want to comment about the relevance of 
our calculations to Navier-Stokes turbulence. If we sub- 
stitute blindly ^ = 4/3 in our results, we predict the 
exponents 2/3, 1.252 26, 2.019 22, 4.048 43, 6.068 60 
and 8.083 37 for £ = 0,2,4,6,8, and 10, respectively. It 
would be tempting to propose that similar numbers may 
be expected for Navier-Stokes and indeed for i = and 
2 this is not too far from the truth. We cannot, how- 
ever, state with confidence that the genuine nonlinearity 
of Navier-Stokes does not change these numbers signif- 
icantly. More work is needed before we can draw final 
conclusions on the rate of decay of the high sectors of 
anisotropy in Navier-Stokes turbulence. 
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APPENDIX A: GAUSSIAN INTEGRATION BY 
PARTS 



The field w{x,t) as well as the forcing are Gaussian 
white noises. This enables us to express T"^{r) and the 



correlation of the force in terms of C°'^{r) and F°'^{r). 
One way to accomplish this is by using the Gaussian in- 
tegration by parts method Q . Using the basic formula 
for Gaussian integration by parts, we get for the third 
moment 



{v°'{x + r,t)wi'{x,t)v'^{x,t)) = 
Jdt' Jdy {w^\x,t)w''{y,t')) 
Sv'^ix + r^t) p 



(Al) 



Sw''{y,t') 



+ {v°'{x + r,t) 



i-^v^\x,t) 



6v'^{x,t) 
5w''{y,t') 



To find out the functional derivative, we formally inte- 
grate v°'{x, t) 

v°'{x,t)= [ dt'dt'v"{x,t') 



dt' w{x,t'Yd^,v'^{x,t') 

dt'[dy[d''drG{x~y)]w^{y,t')d^v^{y,t') 
+ [terms that are independent of w] , (A2) 

and thus 

^^1^ - e{t y)d.v-{y, t') (A3) 

-f [d^drG{x-y)]d,v-{yX) . 



When we plug this result back to Eq. (Al) we face the 
problem of evaluating the step function 9{t — t') a.t t ~ t' 
due to the delta correlation in time of w{r,t). To solve 
this problem in a "physical" way iQ, we approximate 
the delta function of the white noise with a sharp even 
function, perform the integral, and only then take the 
white noise limit. Doing so we obtain the formal result 
9{0) = 1/2 stemming from the fact that we approximate 
a delta function with an even function. Finally, we re- 
mark that this derivation corresponds to the Stratonovich 
interpretation of the stochasic equation Eq. (2.1). 
Next, we perform the spatial integration, arriving at 

(v^ix + r, t)w^'{x, t)v'^{x, t)) 

X jdyG{r-y)dr[K^^%y)d,C^^{y) 

^ ^dyd^drG{y)\K^^^y)d\y^C^^{r-y) , (A4) 



and therefore 
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T"^{r) = dj^^ {v^ix + r)w''{x)v'^{x)) (A5) 
+ iaf,) [dyG{r~y)dr[K>^''{y)d^d.C^^{y 



- jdyd^drGiy) K^^%y)d^y^d^y^C''^{r-y) 



APPENDIX B: SIMPLIFICATION OF THE 
"NONLOCAL" TERM 

Fo r ^ ^ 1, the nontrivial integral on the LHS of 
Eq. ( |2^ ) can be further simphfied. To see that, let us 
denote it by Int and rewrite it (omitting the —1/2 factor) 
with the explicit forms of the Kraichnan operator and of 
the Green function 



/„t = -(e + 2)£ jdyU^dr[y^d^C"^{r-y)\ (Bl) 
jdyU^dr[y^-^V^y''d^,d,C^-{r~y)\ . (B2) 



It is easy to verify that the tensors 9^(7""^ (r — y) and 
y'^y'^d^d^C"'^ {r — y) are divergence- free in both indices 
due to the fact that C°''^{r) itself is divergence free. 
Therefore, to simplify the integrals in Eq. (Bl), we con- 
sider the generic expression 



(dy-d^dry^X'^^iy) , 

J y 



(B3) 



where X°''^(y) is some divergence-free tensor and A is 
an arbitrary exponent. We also assume that A is such 
that the integral is convergent and integration by parts 
is allowed. Then we may write 

fdy-df'dry^X^^y) (B4) 

J y 

= Jdy [a(A - 2)j/^- V2/rX«^(y) 
+ Xy^-^X°''^{y) + Xy^-^yrd'^X^'-'iy) . 
The last formula can be simplified by using identity 
= Jdyd^y^-^yrX'^^iy) (B5) 
= (A-3) Jdyy^-'y^yrX'^^y) 

+ jdyy^-^X^^{y) + j dy y^-^y^d^ X^^ {y) , 
from which we get that for A ^ 3 

(dy-d^dry^X'-^iy) (B6) 
J V 

= Jdy [y^-'X'^0{y)+y^-'yrd0X^^y)] . 



For A ^ 1, we can even do better using the identity 
= Jdydry^-^d^X^^y) 

= (A - 1) Jdyy^-^yrd^X'^^y) , (B7) 
which finally brings us to 

jdyU^dry'-X^^iy) = j dy y^-^ X'^'' {y) . 

(B8) 



Let us now apply Eq. (B8) to the integrals in Eq, 
Assuming that ^ 7^ 1, we get 

D^^ + 2 



47r(C - 3) 
Di{^ - 2 



Jdyy^^^d^C'^ir-y) (B9) 
dyy^-^y'^y-'d^d^C^ir-y) . 



- 5) 

To continue, we wish to turn the second integral into 
the same form of the first integral. To accomplish that, 
consider the following identity: 

= Jdyd^'[y^-^y''d^d,C''''{r-y)] (BIO) 

= (e-3) ldyy^-'y^y-'d^d,C"''{r-y) 

+ Jdyy^-'d'C"^ir-y) 

+ jdyyi-^y-d^d^C^^ir-y) , 
which gives us 

ydy2/«-V2/"9^a,C"^(r-y) 

= jdy [y^-^d^C"^{r-y) 

+ y«-3y''a^a2C"'3(r - y)] . 
Additionally, we have 

= Jdyd.yi-Vd'C'^^iy) 

= e Jdyyi-^d^C''^{y)+ Jdy y^-^y'^d.d^C^f'iy) , 
and so finally we obtain 

J dy y^-'y^y-'d^d^C'^r - y) (B13) 



(Bll) 



(B12) 



Jdyy^-^d^C^ir -y) 



Substituting this into Eq. ( |B9| ), we arrive at the final 
result for Int- 



Int — 



(e-3)(e-5) 



JdyG{y)yi-'d'C^^{r-y) 



(B14) 
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APPENDIX C: CALCULATION OF ^(A;^,^) 

The prefactor A{X;£,£^) was defined by 

A(A;^,C)r^+«r,„,(f) (CI) 
= JdyG{r -y)\r- y\^-^y^Y,^(v) . 

Due to the isotropy of the integral, it is m indepen- 
dent and therefore we speciahze to the m = where 
the Yim{y) is proportional to the Legendre polynomial 
Pe{y ■ z). Setting r = z, the unit vector in the z direc- 
tion, we write the integral in spherical variables (y, 0, 0) 
and perform the trivial (j) integration (for r = z, the 
integrand is independent of </>). We arrive at 



1 



dyy^+^ I d{cose) 



(C2) 



X (y"^ - 2ycos6' -f l)~Pf(cos( 



Using the standrad tricks of Feynmann integrals, one 
can express this integral (at least in the £ — case) in 
terms of gamma functions. Here, however, we choose 
to calculate this integral directly by a strightforward ex- 
pansion of the integrand. This procedure underlines the 
connection between the pole structure and the anisotropy 
label e. 

Let us therefore turn the annoying (y^ — 2y cos 9 + 
l)(?-3)/2 iQYm into a Taylor series in 2ycos9/{l + y"^), 

1 ,22^ (-1)"2" 



2 ^ n! 

n=0 



2 J \ 2 



(C3) 



where B(X,n,^) and C{n,£) are two one-dimensional in- 
tegrals that are given by 

B{n,£) = J rf(cos6l) cos"6'P<?(cos6') , (C4) 

C{X,n,0^ dyy'+''+\l+y')'^-" . (C5) 
Jo 

Before calculating B{n,£), we notice that it vanishes for 
n < £ and for ns that are of different parity than £. The 
first observation is a simple manifestation of the orthogo- 
nality of the Legendre polynomials intimately connected 
to the orthogonality of different irreducible representa- 
tions of the S0(3) group. For n > £ with the same parity, 
we can use the well-known identity 



(C6) 



from which we get, after simple integration by parts, that 
B{n, i) is given by 



B{n,£} ^ (-1) ^j-^i (C7) 



dxx"'-'^{x^ - lYdx 



The above integral can be done explicitly leading us to 
the final result 



B{n,£) = 



i{n-l)---{n-e+l) 
¥11 
£\ 2(-l)'= 
k) n + 2k-£ ' 



(C8) 



E 

k=0 



We now calculate the second int egr al by dividing the 
integration regime [0, +oo] in Eq. ( |C5| ) into a < ?/ < 1 
part and a 1 < y < cxd part. In each part we expand the 
(1 -I- y2')(4-3)/2-n y g^j^fj i^y^ r csp cct Ivcty, aud 

perform the integration. After adding up these terms 
again, we arrive at the following sum 



1 /e-3 



fc! 

fc=0 

X I -2 n - 1 



(C9) 



e-3 



n — k + I 



X + 3 + n + 2k X + ^-n-2k 



Plugging these results back in Eq. (p3D, we get 



OC OQ 



A{X;£,0^- 



2 ^ ^ n\k\ 

n=i,l+2,... k=Q 



(CIO) 



e-3 



1 



n - k + 1 
1 



X + 'i + n + 2k X + i-n~2k 



B{n,£) 



To expose the structure of the poles in the above for- 
mula, it is convenient to change the order of summation 
by defining the new indices q,j, 



n - £ 

q = ^^+k = 0,1,2,, 



2 

n-£ 



(Cll) 



= 0,1,2,. 



and obtain 



(-^V °° 
A(A;£,0 = -^E«^(^'^) 



(C12) 



q=0 



1 



X + 3 + £ + 2q X + ^-£-2q 



where aq{£,£,) are given by 
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3=0 



ii + 2jy.iq-jy. 



(C13) 



^-3 



e-j-q+l Bi£ + 2j,£) 



Equation C12 shows that A{\:i,S^) is given as an in- 
finite series of poles in A. For ^ > this series can be 
shown to converge although for small values of ^ the con- 
vergence is very slow. In the special case of ^ — 2, the 
series is truncated after the first pole. To see why this is 
so, return to the original definition of A{X; i, ^) and set 

A{X;e,2y+^YeUr) = j dy G{r - y)y^Y,^{y) . (CM) 
But since 



v^Yt„,{y) 



1 



{\ + 2-l){\ + i + l) 

(C15) 

then from the definition of the Green function we get 

1 



A{\-£,2) 



(A + 2-^)(A + 3 + ^) 



(C16) 
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